'Tail-end' Bondi-Hoyle Accretion in Young Star Clusters: Implications for 

Disks, Planets, and Stars 

Submitted to AJ, 7-Nov-2007 

Revised, l-Apr-2008 

Accepted, 2-Apr-2008 

bo 

Henry B. Throop 

5_^ Southwest Research Institute 

-^ Department of Space Studies 



Or 



% 



1050 Walnut St, Ste 300, Boulder, CO 80302 
throopOboulder . swri . edu 



John Bally 

C3 \ Center for Astrophysics and Space Astronomy 

University of Colorado 

> 

00 ! UCB 389, Boulder, CO 80309-0389 

m 

o 

Tf 

o. 

Young stars orbiting in the gravitational potential well of forming star clusters pass 
through the cluster's dense molecular gas and can experience Bondi-Hoyle accretion 



ABSTRACT 



from reservoirs outside their individual protostellar cloud cores. Accretion can occur 
for several million years after the stars form, but before the cluster disperses. This 
accretion is predominantly onto the disk and not the star. 

N-body simulations of stars orbiting in three young model clusters containing 30, 
300, and 3000 stars are presented. The simulations include the gravitational potential 
of the molecular gas which smoothly disperses over time. The clusters have a star 
formation efficiency of 33% and a radius of 0.22 pc. We find that the disks surrounding 
solar-mass stars in the N=30 cluster accretes ~ 0.01 Mq (1 minimum-mass solar nebula, 
MMSN) per Myr, with a la width of 50 times due to variations in initial stellar positions 
and velocities within the cluster. The accretion rate scales as M 2 ' 1±01 for stars of mass 
M. The accretion rate is ~ 5 times lower for N=3000 cluster, due to its higher stellar 
velocities and higher temperature. The Bondi-Hoyle accretion rates onto the disks 
are several times low er than accretion rates observed directly onto young stars (e.g., 



Muzerolle et al.ll2005l ): these two accretion rates follow the same M 2 behavior and may 



be related. The accreted disk mass is large enough that it may have a substantial 
and unappreciated effect on disk structure and the formation of planetary systems. 
We discuss a variety of implications of this process, including its effect on metallicity 
differences between cluster stars, compositional differences between a star and its disk, 
the formation of terrestrial and gas-giant planets, and isotopic anomalies observed in 
our Solar System. 



Introduction 



Stars form from the collapse of dense cores in giant molecular clouds (GMCs). While some stars 
form in relative isolation or in small groups, the majority of stars in the Solar neighbo rhood appear 



to be born in transient c lusters containing hundreds or thousands of members (e.g., lAdams et al 



20061 ; lLada fe Ladall2003l ) . Dense cores have masses ranging from under IMq to over 10 Mq for the 
most massive cluster-forming environments. Self-gravity and efficient cooling by dust and molecular 
line radiation leads to collapse and fragmentation. The timescale for the birth of an isolated star 
can be estimated by considering the accretion rate onto a protostar from an isothermal sphere: 
Tacc ~ M/M ~ GM/(? S , where M is the final stellar mass and c s is the effective sound speed in 



the core. Assuming a 1 M & star in a cloud with c s = 0.35 km s , we find that r a 



10 5 yr. 



In contrast, the observed age spread of young stellar objects (YSOs) in clusters indicates that the 
formation time-scale for an entire cluster is a few million years. Thus, only about 3 to 10% of 
the final population of YSOs in a cluster are expected to be in their main accretion phase (Class 
/ Class I) at any time. Star formation then can be characterized by a hierarchy of at least 
three timescales: the 10 5 yr scale to form an individual star; the 10 6 yr to form a cluster; and the 
10 7 yr timescale of an OB association. Planet formation operates on roughly comparable timescales: 
10 5 yr to form meter- to km-sized bodies; 10 6 yr to form planetesimals and cores (for the terrestrial 



and giant planets, respectively); and 10 7 yr to complete the planets (ILissauer &: Stevenson 



Nagasawa et al.l 120071 : 



Lada fc Ladall2003l : lThroop et al.ll200ll : lBallvll200ll : IWeidenschillind 119971 ) . 



2007 



As the cluster gives birth to more stars, the oldest stars will evolve into Class II and III YSOs 
surrounded by protoplanetary disks. For at least a few million years, as these YSOs move through 
the cluster, they have a chance to accrete additional material from reservoirs of dense gas remaining 
in the re gion. Bondi-Hoyle (BH) accretion describes the rate at which material will be added to 



the star (Bondi 



1952 



Bondi fe Hovldll944l ). The accretion rate is generally much less than that 



experienced b y stars in the collapse and accretion phases of their formation: 10 8 -M© yr 1 vs. 



10~ 5 Mq yr -1 ( Bate fc Bonnellll2005l ). The total mass accreted is small compared to the star, but 
it can be large compared to the mass of the protoplanetary disk. Moreover, the disk intercepts 
material falling toward the star, so mass is deposited onto the disk and not the star. Thus, this 
'tail-end accretion' can have profound consequences for the evolution of very young planetary and 
pre-planetary systems. Because accretion scales as M ex M 2 , it is likely to be most important for 
the larger members of a forming cluster that spend the longest time moving through dense gas. 



Typical star-formation efficiency (SF E) in a cluster-forming cloud core is 10% — 30% (e.g., 
Jorgensen et al.l 120071 ; lLada fe Ladall2003l ). The majority of the gas is not consumed by stellar 
birth, but remains in the cloud for millions of years, even in the presence of ongoing star formation. 
For instance, dense molecular cores located as close as 0.1 pc to the ONC cluster core continue 
to produce stars; the OMC-1S core 90" southwest of the Trapezium stars contains several dozen 
currently accreting protostars (jZapata et al.ll2006l ). Likewise, the moderate mass cluster IC 348 
in the Pers eus Molecular cloud contains more than 400 young stars with an age spread of least 
2.5-5 Myr (JMuench et al.l 120071 ) - yet, this cluster contains at least 20 cla ss /I protostars , many 
of which drive outflows at it southern periphery (jWalawender et al.l 120061 ) . iMuench et al.l ()2007l ) 
propose that this cluster was built over several crossing-times and that the current generation of 
class 0/1 YSOs will eventually fall into the cluster core and merge with it. 

Stellar winds, supernovae, an d other processes ev entually remove the gas; clusters with gas re- 
maining beyond 10 Mya are rare (jLada fc Ladall2003l ). When only low- to intermediate-mass stars 
are formed, proto stellar winds and jets are the most likely agents responsible for the removal of re- 



maining gas (e.g. JWalawender et al.ll2005l ). In massive cloud cores, outflows may not inject enough 
energy and momentum to disperse the gas, enabling star formation to continue. Furthermore, 
winds may increase the turbulent motions and inter nal pressures of the remaining gas, thereby 
increasing the tendency to form more massive stars (|Tan fe McKed 120041 ) . When massive O and 
B stars form, their intense UV radiation fields terminate star formation by heating and removing 
the remaining gas. UV heating and ionization will bring star formation to a halt in a crossing time 
which can be estimated from the radius of the region, L, and the sound speed in photo-ionized 
plasma, c s ~ 11 km s . Typically, this time-scale ranges from 10 5 to a few times 10 6 years. In 
some cases, UV radiation may also trigger star birth: as expanding ionization fronts overrun and 
surround cores, the rocket effect of photo-ablating plasma combined with the high pressure of the 
H II region can result in radiation-driven implosion, gravitational collapse, and triggered star birth. 

Star formation is observed over a broad range of size scales. The Orion nebula is the proto- 
typical example of a nearby massive region that h as spawned over 2000 stars in the last few Myr 
([Hillenbrand &: Carpenteril2000l : lHillenbrandlll997l ). The dynamics of the gas is now dominated by 
UV radiation and winds emerging from the massive Trapezium stars at its core, but for several Myr, 
the region was cool and dark as low-mass stars formed and orbited through the cluster. The Orion 
nebula's solar-mass stars are observed to h ave ages of up to 1-2 Myr, far older than the 10 4 -10 5 
year- old Trapezium stars (Bally et al. 19981) . The NGC 1333 region in the Perseus molecular cloud 
(e.g., I Jorgensen et al.l 120071 ; lLada et al.l 119961 ) provides a good example of a medium-sized cluster 
that has formed about 150 low- to intermediate-mass stars within the last 1-2 Myr. Protostel- 
lar winds appear to dominate the gas dynamics and may regulate the rate of star formation. At 
small er size scales, star fo rmation in Taurus is divided into numerous small, distinct clumps of ~ 15 
stars (JGomez et al.l ll993). These small clusters have far lower spatial density than Orion, but with 
similar SFE. 



Across the range of stellar number, the timescales and spatial scales of star formation are 



remarkably uniform: all three have ages (and lifetimes) of several Myr, and diameters 0.5 — 1 pc. 
Large regions are characterized by higher gas and stellar densities, but only slight increases in 
spatial scale: for instance, the width of the Orion nebula c ore is ~0. 1-0.2 pc, comparable to the 



low-m ass cores in Taurus with their 100-times lower density (jjohnstone &: Ballvlll999l : lOnishi et al 
19961 ). 



This paper reports on numerical simulations of Bondi-Hoyle (BH) accretion onto stars and 
disks during the several Myr after stars have formed, but before clusters have dispersed. This is 
assumed to be a distinct, long-lived phase after the initial mass function (IMF) has largely been 
established. Accretion onto the star /disk continues, but at a rate far lower than any accretion 
during stellar formation. This 'tail-end' phase of accretion has received little attention, but may 
prove to have substantial consequences, because it is during this 1-10 Myr period when the disk 
exists that planetesimals and young planets complete all but the final stages of their formation. 

BH accretion is a well-studied process with several applications to stellar formation and evo- 
lution. It has been studied for its role in determining the stellar IMF. The origin of the IMF is not 
well understood, with debate primarily between two models. In the first, stars fragment from the 
cloud with only a small seed mass, and acc rete the rest of their material via BH 'competitive accre- 
tion' from the cluster over a few 10 5 yr (cf. iBonnell et al.ll2007l . and referenc es therein). This mode l 
is supported by extensive numerical simulations; for instance, the work of iBate fc Bonnelll (120051 ) 
used SPH calculations to model the collapse of a cloud with radius of 0.19 pc and mass 50 M Q , for 
~ 300 kyr, during which the IMF was almost completely determined. The second model proposes 
that the IMF is set at the time of cloud fragmentation. A molecular cloud breaks into a spectrum 
of pr e-stellar cores, and the se cores condense into stars without accreting substantially more mass 
(cf. iKrumhorz et alj|2005bl . and references therein). This model is supported by direct millimeter 
observations of young cores, which show pre-stellar con densations in a distribution that clo sely 
matches a Kroupa-type IMF over the range 0.05 - 5 M ? ()Andre et all 120071 : iMotte et al.lll998l ). 



Recently, iPadoan et al.l (|2005l ) invoked BH accretion to explain the observed accretion signa- 
tures onto young stars, which we address later in this paper. On longer timescales of 10 10 yr, BH 
accretion within gl obular clusters has been used to explain the 'self-pollution' o f high-metallicity 



veneers onto stars (Thoul et al. 



2002; 



Parmentier et al. 



2000 



Jehin et al.l ll999). As a historical 



note, an early paper by iLvttletonl (|196ll ) proposed that the entire Solar System was accreted by a 
BH process as the Sun tr aveled through thousand s of cold molecular clouds, but such a scenario 
was shown implausible bv lAust fc Woolfsonl (|1973l ). 



In this paper, we use N-body simulations to model stellar dynamics in three clusters containing 
30, 300, and 3000 stars, and containing a large reservoir of gas which disperses on timescales of 
several Myr. This paper focuses on the roughly solar-mass stars since they are the ones of greatest 
interest to studies of planet formation. The accretion of gas onto the stars and their disks is 
calculated using the Bondi-Hoyle equations, with no new hydrodynamic simulations. In § [2] we 
give background and motivation. Numerical simulations are detailed in § [3l and § 2] discusses the 



results. We compare to observations in § [5], and we discuss implications for the formation of stars, 
disks, and planets in § [6J 



2. Background and Motivation 
2.1. Bondi-Hoyle accretion 



The Bondi-Hoyle accretion (BHA) rate can be written as (e.g.. lEdgarl 12004 ) : 

M^ = Txn^m H vR\ cc , (1) 

where n is the gas number density, [xmu is the mean molecular weight, v is the stellar speed 
relative to the gas, and Racc is the accretion radius, which is the distance outward to which the 
star's influence on the gas can be felt. For an isolated star, i?Acc is the BH radius Rb given by 

2GM 

<y + cj) 

where G is the gravitational constant, M is the stellar mass, and c s is the sound speed. For stars 
in a dense cluster, accretion can happen only out to the tidal radius Rt, where 

a ( M \ 1/3 

for a star at distance a from the cluster core, enclosing a cluster mass M enc within its orbital 
dista nce. The star accret es material out to Racc, which is the lesser of the gravitational and tidal 



radii (jBonnell et al.ll200ll ): 



R Acc = mm(R Bl R T ) . (4) 

For the parameter space explored here, Rb <C Rt in almost all cases. In our simulations, the two 
become comparable only for short periods at the apoapses of highly elongated orbits. 

For a typical cool molecular cloud of 25 K and molecular weight 2.3 run, the sound speed is 
c s ~ 0.3 km s . Using n = 10 4 cm -3 and v = 2 km s _1 , we find Rb ~ 500 AU, and the accretion 
rate is Mb ~ 10 -8 M & yr _1 . For most YSOs in a cluster, Rb will be larger than typical disk sizes 
of 10-200 AU, so the gravitational cross-section for accretion is many times higher than the YSO's 
physical cross-section. 

The velocity v is the stellar speed relative to the gas. v is difficult to calculate directly because 
the gas is turbulent and characterized by many velocities on different scales. A simplification could 
be made that the gas is at rest, but this would result in a velocity estimate that is usually too 
low, causing an overestimate of the accretion radius and the accretion rate. We use two different 



methods to estimate v. Both are based on the assumption that the gas can be thought of as being 
composed of small clumps that travel primarily on gravitational orbits through the cluster region, 
similar to the stellar motion. This assumption is supported by observation s of the motions o f 



pre-stellar condensations in young clusters, moving at roughly virial velocities (| Andre et al.l 120071 ). 
Therefore, in our first approach, we estimate v based on the stellar velocity v* and virial velocity 

v v by 

v = {vl + vW (5) 

where 

v v = [2 (KE + PE)/M] 1/2 . (6) 

KE and PE are the total kinetic and potential energy of the system, and M is the total mass. The 
virial velocity evolves with time, and generally gets smaller as the system evolves. 

The second method also assumes that the gas clumps behave dynamically like stars, and thus 
calculates v by recording the actual velocity v n between a star and its nearest stellar neighbor, at 
each timestep, and then taking 

v = v n . (7) 

In this way the distribution of velocities between stars and gas clumps orbiting within the cluster's 
potential well can be estimated. A byproduct of this method is that transitions between velocities 
are not as smooth as they would be in reality; that is, when a star's nearest-neighbor changes, its 
v n also immediately changes. For conceptual simplicity some of our plots therefore use eq. [5], but 
all of our final results use eq. [7] . 

Bondi-Hoyle accretion is a simple and robust process, but it can be affected by several environ- 
mental factors. St ellar winds from Class II YSOs or T Tauri stars have M w ~ 10~ 8 Mq yr" 1 and 
v w ~ 200 km s _1 ( Rav et al.l 120071 ). The outward ram pressure from such a wind (P = Mv/R 2 ) 



could be enough to prevent accretion from a competing inward BH flow, because the mass rates are 
comparable, and the wind velocities greatly exceed the BH velocities. However, these winds and 
jets tend to be quite narrow, and can only stop accretio n for material that passes directly through 



their paths. HST observations by lHartigan et al.l (|2004l ) of several young T Tauri stars found the 



jets to have opening angles of ~ 20° at distances 20-50 AU. At greater distances, the opening angles 
are smaller and the jets become essentially collimated. Thus, near the stars the winds cover solid 
angle Q < 1 sr, and this decreases substantially at Rb- Stellar winds may therefore be responsible 
for reducing BH accretion by a fraction $7/(4-71") ~ 0.1. Conceivably, if the flow passes through the 
wind cone several times this fraction could be increased; however, we argue below that the BH flow 
accretes onto the disk rapidly, and once on the disk it is certainly immune from polar outflows. 
The detailed nature of the outflow- inflow interaction could be studied with future simulations. 



Edgar fc Clarkd (|2004j) showed that radiation pressure from low-mass stars is too weak to affect 
the accretion rate, but that above IOMq, radiation pressure becomes strong enough to inhibit BH 
accretion entirely. The transition between these two regimes is sharp. Because we focus on solar- 
mass stars and our models have very few high-mass stars, this effect has no consequence on our 
results. 

BHA is normally modeled as accretion onto a point source or star. The flow is gravitationally 
focused from the initial width of Rb to a much smaller size as it nears the star (Figure [1]). The 
flow does not accrete onto the star's front hemisphere, but continues beyond the star, where it loses 
much of its energy when it collides with opposing streamlines. The flow behavior can be modeled 
analytically for subsonic flows, and has been th e subject of numerous 2D and 3D s i mulations for 
higher-speed, turbulent f lows at Mach 3-10 (e.g., 



Foglizzo fc Ruffert 



1997 



Ruffert 



1997 



Krumholz et al. 



2005a 



Edgar 



2004 



Ruffert 



1999 



Benensohn et al.lll997l ). These simulations find that as the 



accreted material collects behind the star, it sets up a transient, tenuous 'accretion disk' of size 
~ 0.2 Rb- The disk oscillates unstably between prograde and retrograde rotation (the so-called 
'flip-flop' instability), and eventually accretes onto the star on Keplerian timescales. The accretion 
flow has rate variations of a factor of 2-5, and the tenuous disk has mass ~ 10 -3 MMSN, for 
conditions in our model cluster. Most of the mass contained inside Rb is ultimately accreted onto 
the star. Angular momentum is also accreted onto the disk, but less efficiently: a gradient of a few 
percent i n the backgro und gas velocity can inhibit most of the angular momentum accretion onto 
the star (JRuffertlll999h . 



These simulations have been performed using stars as accretors, but to our knowledge there 
are no published numerical simulations of BHA onto stars with pre-existing disks. The existing 
work shows that the accretion flow is not directly onto the star, but passes through the plane 
of the disk many times before being accret ed onto the star. However, upwards of 80% of young 
stars in clusters are observed to have disks (ISmith et al.ll2005l ). and for the MMSN case, the disk 
density is on the order of 10 3 times higher than the accretion flow density. Infalling material will 
be rapidly intercepted and absorbed by the disk. Only in the relatively unlikely circumstance that 
the star moves through the medium with its disk aligned precisely edge-on, or with no disk at all, 
can accretion directly onto the star be significant. 

In some ways the accretion flow onto th e disk can be compared to the proto-Jovian nebula 
'starved disk' model of ICanup fe Ward! (J2002J), where gas accretes from the solar nebula disk onto 
the circum-Jovian disk orbiting within it. The two disks are co-planar and move at constant velocity 
to each other, making mass transfer between them a smooth, linear process. The molecular cloud 
case is much more complex, because the orientation and velocity of the disk relative to the ISM 
is constantly changing. The accretion flow will almost always have a different angular momentum 
vector than the star's disk. In the short term, accretion onto the disk will always result in mass 
gain. However, if the incoming angular momentum is opposite that of the disk, one can imagine 
that orbital decay and mass loss onto the star will soon follow as the disk is perturbed. There are 
thus two distinct accretion rates: one from the ISM onto the disk, and one from the disk onto the 



star. The two are related, but the disk acts as a buffer between them, on the viscous timescale of 
10 5 — 10 7 yr. Our current understanding of the interaction between the accretion flow and the disk 
is very limited; we discuss the relationship between the two accretion rates further in § 15.11 

The numerical simulations presented here include large-s cale density gradients in the back- 



ground gas field, but do not include gas turbulence. Work by iKrumholz et al.l ([2006) generalized 
BHA to a supersonically turbulent medium. At high turbulence, the accretion rate is limited by 
large-scale vorticity, which can suppress accretion. In general, they found vorticity to be important 
in the cases Rb/1 > 1, where I is the largest size scale of the turbulence, comparable to the distance 
between protostellar cores. In the core of our densest simulations, ra* = 10 6 pc -3 , so Rb/1 ~ 1/4. 
This can cause a moderate reduction in accretion; however, such dense regions occur only in our 
largest simulation, only at its core, and only at the beginn ing. In the vast majori ty of our parame- 



ter space, Rb/1 <C 1, and the BH equations remain valid. IKrumholz et al.1 (|2006l ) argued for some 
inhibited accretion in molecular clouds, but they assumed far larger Rb than used here (20,000 AU 
vs. 500 AU), due to their initial conditions of colder gas and lower stellar velocities. 



3. Numerical simulations 

We use N-body simulations of three model clusters to trace the dynamical and environmental 
histories of individual stars. Each cluster contains a background distribution of gas from which 
these stars can accrete. The code NBODY6, a sophisticated Hermite-method integrator designed 
for accuracy in handling of close approaches between stars, is used to simulate the evolution of each 



star cluster (lAarsethlll999l ). NBODY6 is fast and proven with a long heritage. It allows for flexible 
control over initial input conditions, and enables the inclusion of a background molecular cloud 
whose gravity may dominate the motion of stars. The cloud is modeled as a gravitational Plummer 
potential of a given width and mass, which decays smoothly after an initial delay. NBODY6 treats 
the gravitational effect of the gas on the stars, but the stars have no direct feedback on the gas: 
the cloud shape, density structure, and evolution are directly controlled only by input parameters. 

Our calculations are performed in two discrete steps. First, N-body simulations of a cluster 
generate the positions and velocities of each star, outputting at fixed intervals. Second, these 
quantities are used to calculate Mb at each step for each star, using eqs. [IHZ]- The output timestep 
used is 2 x 10 4 yr, which is sufficient for convergence. 

Several simplifying assumptions have been made. First, the simulations are not entirely self- 
consistent in that mass accreted onto stars is removed only as a general global process, and not 
locally in response to accretion. However, mass lost to accretion by the end of the simulations is 
less than 1% of the original cloud mass, so ignoring this has only a minor effect on stellar motion. 



Rufierti ( 1996 ) shows that the drag on an accreting body is approximately M v, so the stars' velocity 



change is also insignificantly small. Because the sound speed is comparable to the stellar velocities, 
any 'tunnels' created in the cloud by BH accretion will be rapidly filled, making a local treatment 



unnecessary. Second, the increase in mass of individual stars due to BH accretion is ignored during 
the N-body simulations; it is only calculated afterward. The increase in stellar mass is usually 
< 1%, making only small changes to a star's motion. This assumption may underestimate growth 
for some high-mass stars, by allowing them to retain their full velocity even as their mass increases. 
However, our emphasis is on the consequences of BHA for solar-mass stars, and ignoring this small 
effect sets a conservative lower-bound. 

For simplicity, it is assumed that all stars are formed at t = 0. Observations show that star 
formation is an on-g oing process that beg i ns as the cloud collapses and ope rates continually until 
the cluster disperses dSwift & Welch! boOTJ : lAUen et alJbo07l : lTan et alJbood . e.g.,). Therefore, our 
model overestimates the total mass accreted by the latest-formed stars, by assuming that all stars 
form at the start. However, the rate of accretion onto a star is not affected by this simplification. 
We implicitly assume that the st ar formation timescale is much shorter than the cluster lifetime 



(El megreenll2000l ; lHartmann 



2001 



but also see Tan et al. 



2006 



Tassis fc Mouschoviasl 12004 ). Stars 



are formed with a pre-set IMF; any mass evolution of the IMF due to ongoing accretion naturally 
shows up in our results. 

We take the initial stellar and gas distributions to be three-dimensional Plummer spheres, with 
a profile given by 



p s (r) oc [1 + (r/r ) 2 ] 



-5/2 



(8) 



where r is the distance from the cluster core, and r$ is a measure of the cluster's size, <±v7o 
of the mass is contained within this distance. This symmetrical distribution idealizes the shape 
and radial distribution of real cloud cores. Altho ugh some regions (such as Orion) are highly 
asymmetrical, observations by IShirlev et al.l ( 20031 ) found a majority of dense cloud cores (10 2 - 
10 4 M & ) to be symmetric, with aspect ratios < 1.3. For cores such as these, Plummer spheres 
provide a remarkably good fit to the radial distribution over a variety of size ranges. Figure [2] shows 



the lin e-of- sight density meas ured in four optically thin condensations measured by lBeuther et al 



(J2002I ) and I Andre et al.1 (120001 ) . The largest three of these are young high-mass star- forming regions 
in Cassiopeia, with masses 1200-2300 Mq, while the smallest is an individual pre-stellar core in 
L1689. All four have flat inner regions with steep outer edges, and all are readily fit by Plummer 
profiles (eq. [U integrated line-of-sight). Because of the broad appropriateness of the Plummer 
distribution, and its easy analytic form, we have adopted it for our initial distributions. An improved 
future model could explore the effects of an asymmetrical cloud core shape. 

The Plummer sphere is a smooth function that does not have the local density variations of 
a realistic molecular cloud. Its density changes only with radial distance. However, because the 
accretion rate for an individual star depends linearly on density (eq. [1]) , local variations in density 
will in general average themselves out over an orbital period. 

Our simulations do not consider the effects of binaries. Long-lived binaries are almost always 
the result of primordial binaries that have been hardened due to subsequent three-body interactions 



10 



(e.g.. iGoodwin et al.ll2007l ). We include no primordial binaries, because we are mostly interested 
in the evolution of single stars around which planetary formation is best understood. For binary 
stars separated by much less than Rb ~ 500 AU, accretion will depend only on the total mass of 
the system, regardless of whether it is a single or multiple star. 

A s with previous simula tions, we ignore magnetic fields. The effect on accretion is most likely 
small (JKrumholz et al.ll2006l ). but requires further study. 



3.1. Initial parameters 

The parameters for our three simulations are listed in Table [TJ The small cluster model is a 
small, dark cloud, similar to one in Ta urus, with N = 30 stars and a gas mass of M g = 30 Mq 



0- 

Masses are distributed using the IMF of lKroupa et al.l (119931 ). The minimum stellar mass is 0.1 M®, 
and the mean is 0.5 M®, giving a total stellar mass of M s = 15 M®. The maximum is determined 
by the IMF, and is 1.2 M®. Star formation efficiency (SFE) is the ratio M S /(M S + M g ) and is 
33%. The gas is retained entirely until t = 2 Myr, after which it is smoothly lost to simulate the 
dispersal by stellar winds and other processes. The simulation ends at t = 4 Myr, which is picked to 



corre spond to typical estimates for the lifetime of star formation (JBriceno et al.l 120071 ; I Allen et al. 
20071 ). 



The medium and LARGE models assume the same SFE, gas radius, cluster radius, IMF, and gas 
timescales as the small cluster. Parameters that change are the number of stars, the maximum 
stellar mass, and the total system mass. The medium cluster models a region like IC348, with 
N = 500 stars, while the large cluster models the environment in a young version of the Orion 
Nebula Cluster (ONC) with N = 3000, before O/B stars have heated the gas and stopped BH 
accretion. 

The Plummer radius in all cases is taken to be r$ = 0.22 pc, which implies a FWHM in 
column density of 0.3 pc. This size results in properly normalized values for the masses and peak 
densities of clusters: e.g., with this tq, our large case has a peak density of N = 10 6 pc~ 3 and 
n = 10 6 cm" 
models. 



which fit observations of Orion. Peak densities are also correct for the two other 



We take initial stellar velocities to be isotropic, drawn from a distri bution specifie d by g(q) = 
q 2 (l — q 2 ) 7 ' 2 , where q = v/ v P . is the ratio of ste llar to escape velocities (jAarsethl 120031 ) . Stars are 
initially virialized. Work by I Adams et al.l ([200a) explored the effect of initial stellar velocities and 
viriality on cluster evolution. They found that stellar velocities in general became more radial with 
time, especially as gas is removed. We show in £JJ] that most of the BH activity occurs during brief 
passages through the cluster core by stars on radial orbits. Our assumption of isotropic initial 
conditions thus represents a conservative case: choosing subvirial (i.e., radial) velocities would 
increase the number of close passages through the cluster core, increasing BH accretion. 
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The initial radial profi le of the gas envelop e is the same as for the stars. The gas mass then 



decays over time, given by (JKroupa et al.l 12001 



M(t)- t<td : M * (0) N) 

9[) \t>t d : M g (0)/[l + (t-t d )/D] ■ W 

Here t is the time, tj is the delay time before loss begins, and D is the timescale for loss to half the 
initial mass. D is set so that the gas mass at 4 Myr is 5% of its value at 2 Myr. 

We assume two different gas temperature models. In the first, we take T=5 K, T=15 K, and 
T=40 K for the small, medium, and LARGE clusters, respectively. In the second, we use T=25K 
for all cases. Except where otherwise mentioned, all of our results are reported using the first model; 
the second model is used only as a test case to 'unwrap' the effects of cluster size and temperature 
on accretion rate. In all cases, the sound speed c s is given by 

\lim H J 

No new stars are produced during the runs; a few stars are removed as they are ejected beyond 
~ 5 pc from the core. All runs are stopped at 4 Myr; by this time the gas has mostly dispersed, 
and the accretion rate is insignificant. 

We performed multiple runs for each set of initial conditions, varying only by the random 
seed, and combined these together in order to have equal statistical sampling of each cluster. Our 
findings reflect the results from 6000 stars for each case: 200 runs for N=30, 20 runs for N=300, 
and 2 runs for N=3000. 



4. Results 

4.1. Cluster morphology and validation 

Figure [3] shows images of the large cluster at the start and end of the simulation. During the 
first 2 Myr, the cluster expands somewhat, predominantly due to low-mass stars that are given high 
velocities. The expansion speeds up significantly by about 2.5 Myr, when the reduced gas potential 
begins to be felt. By the end, the cluster volume density has decreased by a factor of about 40 
inside of 0.1 pc, and increased slightly outward of 1 pc. About 15% of the stars are unbound 
from the cluster; more are lost after the run has stopped. The SMALL and medium clusters evolve 
generally in a similar way to the large cluster. 

The simulations have been validated for numerical consistency by checking for conservation of 
angular momentum and energy. Energy is conserved up until the point at which gas loss begins. 
After this time, the total energy smoothly decreases to reflect the central potential's reduced mass. 
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Angular momentum is similarly conserved up until gas dispersal, at which point it increases to 
reflect stars' expanding orbital radii. 



4.2. Stellar motions and Bondi-Hoyle accretion 

Figures [4] and [5] show the positions and distances for nine stars: three typical examples from 
each of the three clusters. Each of the stars takes brief, rapid passages through the cluster core at 
< 0.1 pc, but spends most of the time at apoapse near 0.5 pc. For the first 2 Myr, most stars in 
the clusters move on eccentric orbits with e > 0.7, calculated by measuring periapse and apoapse 
on a typical orbit. Stars in the small cluster orbit slowly due to the low central potential, and 
stars have only ~ 1 crossing period before the gas is lost. These stars' orbits are most easily 
perturbed by other stars, because the central gas potential is low relative to the stars' individual 
mass. In comparison, the large cluster has a higher mass and thus shorter orbital periods, so 
stars have ~ 10 crossings during the simulations. These stars follow more regular orbits and are 
less frequently perturbed than stars in the Small cluster. The medium cluster is bounded by the 
other two. Expansion of the clusters is clearly evident in all three runs: when gas begins to be lost 
at 2 Myr, stellar orbits expand, and their periods increase. Because the gravitational potential is 
not a point mass, stellar orbits are not keplerian and tend to precess, most visibly in the large 
cluster. As the gas is lost, orbits become less periodic, and some stars are lost from the cluster. 

Stellar velocities v* for the same nine stars are plotted in Figure [6j Velocities are lowest in the 
small cluster (0.1 — 1 km s _1 ) and highest in the large (1 — 10 km s _1 ), because of the higher 
gravitational potential. Velocities increase at periapse and are slowest at apoapse. Because of the 
dispersed central potential, the velocity change through a star's orbit is not as strong as for regular 
keplerian motion. As the gas disperses, the velocities decrease. 

The accretion rate Mb is also plotted in Figure [6j The two curves of Mb assume two different 
methods (eqs. [3HZ1) of transforming the stellar velocity v* into a stellar-gas velocity v to be used to 
calculate Rb- The curve for v v (virial method) is smooth because the virial velocity changes only 
slowly. On the other hand, v n (nearest-neighbor method) changes rapidly. The values for v n are 
similar to v v , but with a high-frequency noise added due to the rapidly changing velocity of a star's 
actual nearest neighbor. 

Typical accretion rates are 10 -9 to 10 -8 M@ yr _1 in the clusters. The small cluster has the 
lowest gas densities (which decreases Mb), but also has the lowest stellar velocities (which increases 
Mb). These two terms nearly cancel out: accretion in the large cluster is only a factor of 3 lower 
than the small case, even though the gas densities differ by a factor of 100. 

The change in accretion rate through an orbit can be most easily seen in the 'virial' Mb plots. 
The accretion rate for an individual star changes by up to a factor of 10 through its orbit. The 
highest accretion rates usually are near the cluster core, where the density is highest. Velocity is 
also highest at the core, but the effect of the higher density dominates. The accretion rate for all 
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stars drops as gas is lost from the cluster, although the drop in gas density is (once again) partially 
offset by the lower velocity. 

The 'neighbor' M-q curves show the same effect, but the accretion in effect has a noise added to 
it due to the velocities of its rapidly changing nearest neighbors. The accretion peak often occurs 
near periapse, but can occur elsewhere if a star has a fortuitous encounter with a slow- moving 
neighbor star. 

The instantaneous accretion rate Mb can be summed at each timestep to compute the total 
mass accreted for each star. This quantity AMb is plotted in Figure [7] for each of the nine stars. 
By the end of the simulation, the stars have accreted between 0.1% and 9% of a stellar mass. The 
accretion is generally episodic, and is most easily visible in the 'virial' case (grey lines), where the 
velocity is not influenced by the particular nearest neighbor. The episodic accretion is easily visible 
as stairsteps in the accretion. These stairsteps each deliver ~ 0.5-10 Jupiter masses, in ~ 100 kyr- 
300 kyr steps (lMj = 0.001M©). As cluster size increases, the episodes of accretion are shorter, 
but more frequent. There is a diversity in how episodic the accretion is: while the stairsteps are 
visible for some stars, several on more circular orbits show an accretion rate that is nearly constant, 
until it drops at the end. 

For these nine stars, the largest individual mass gains (both in absolute and relative mass) are 
in the SMALL cluster, where one YSO has increased its mass by nearly 9%. Examination of this 
star's particular orbit (Figures [U [5l upper right) indicates that it was formed near the cluster core, 
and stayed inward of 0.2 pc where it was immersed in high-density gas for nearly the entire 4 Myr 
simulation. On the other extreme, one star shown in the LARGE cluster accretes less than 0.1% of 
its mass. This star's initial position is far from the cluster core at 0.4 pc, and it never once passes 
through the center before the gas starts to disperse. 

Figure [8] plots the total AMb as a function of M for every star in the three simulations. Typical 
solar-mass systems accrete ~ 1 MMSN during the runs, with stars in the SMALL cluster accreting 
slightly more than those in the large. A small fraction accrete 10~ 3 M Q or less; these are stars 
that are ejected before even one passage through the cluster core. On the other extreme, about 1% 
of the solar- mass systems accrete as much as 0.1 Mq. The full range of accretion for solar-mass 
stars spans four orders of magnitude! This diversity in accretion amounts may contribute to the 
observed diversity of disk sizes; for instance, in Orion many stars have disks that are not resolvable 
at the resolution limit of ~ 20 AU, while other stars likely of similar mass have disks hundreds of 
AU across. 

A clear trend between AMb and M is visible across a factor of > 100 in stellar mass. Fits 
to the results in Figure [8] are computed using a least-squares method in logarithmic space (i.e., a 
linear fit to the log of the mass and accretion rate) . The mass accreted is listed in Table [2] and is 
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We find that solar mass stars on average accrete slightly less than ~l-4 MMSN worth of 
material, depending on the cluster size. Stars in the small cluster accrete more material because 
of the lower stellar velocities there, which more than cancels out the effect of that cluster's lower 
gas density. Likewise, gas densities in the large cluster are high, but this is more than canceled 
out by the higher encounter speeds. The difference between the large and small clusters is a 
factor of ~ 5, relatively small considering the total range in accretion rate of 10 5 times during the 
simulation. We assume here that stars accrete for 4 Myr; if a cluster is longer- or shorter-lived than 
this, or its accretion is inhibited by the birth of massive stars and an H II region, then the total 
should be adjusted. 

The effect of the different temperature assumptions makes minimal effect on the results (Ta- 
ble [21 bottom). In the small cluster, increasing the temperature from 5 K to 25 K decreases AMb 
by about 30%. In the medium and large models, the difference is <5%; these clusters have higher 
stellar velocities in the first place, so the effect of c s is expected to be much weaker. Therefore, 
most of the difference between the three clusters is a result of the change in density and velocity, 
not temperature. 

Because of the large time- variability of Mb , the mean rate calculated above is not representative 
of the median instantaneous rate that an observer would be likely to measure in a cluster at a given 
instant. Because most of the mass is deposited during brief passages through the dense core, the 
instantaneous Mb should be lower than the mean Mb- We therefore performed a similar fit as 
above, but using the entire ensemble of data points {e.g., 3000 stars at 200 timesteps for 2 runs for 
the large cluster). Best fits to these rates yield 
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These median rates are 2-4 times lower than the mean rates; i.e., an observer would be biased 
to underestimate the net accretion rate, because most of the accretion happens during brief events 
that the observer is less likely to see! 

The plots appear to show that the range in M-q between individual stars at the same mass 
is greater than the systematic variation between the individual clusters. To quantify this range, 
Figure [TOl shows the spread of Mb for just the solar-mass stars in the large cluster. We calculate 
the 1<7 width in median M-q to be a factor of 50; that is, 67% of the instantaneous Mb values for 
solar-mass stars lie within the range 5 x 10~ 10 to 2.5 x 10~ 8 M yr _1 . (Or, at any given time, 
l/6th of the stars in the cluster are accreting more than 2.5 MMSN per Myr!) This range reflects 
changes in accretion rate along a star's orbit, in addition to each star's varying initial conditions. 
The spread is comparable in all three clusters. Thus, wide variations in accretion rate are expected 
to be observed in clusters, even for stars of the same mass, observed at the same time, and even in 
our highly idealized case of a smooth gas distribution. Local variations in gas density (which we 
do not model) would add an additional spread on top of this. 

Figure [11] shows the range in total accretion amount AMb for solar- mass stars, where we 
calculate the la width to be a factor of 7 (i.e., l/6th of the disks accrete more than 4 MMSN during 
the simulation). This range is less than the la width of 50 in Mb because it is time-averaged, but 
still exceeds the difference between the three clusters. 

In all cases, the exponents are very close to the value of 2 for BHA through a uniform medium 
at a fixed speed. This is not entirely surprising, because the orbit of an individual star is virtually 
independent of its mass. That the value is slightly above 2 may be a result of dynamical friction, 
which causes the most massive stars to slow and begin settling toward the cluster center. This mass 
segregation is a well-studied effect on long timescales (JBinnev fe Tremaindl2008l ). A more detailed 
model which included the dynamical effect of increasing stellar mass on the N-body simulations 
(which we ignore) would probably increase this exponent. 



Observational effects 



5.1. Comparison with observed accretion rates 

Accretion directly onto young stars can be detected through signatures in their optical contin- 
uum veiling or Ha emission line profiles. Recent observations have found the unexpected result that 
stellar accretion scales with stellar mass as M oc M 2±a2 , with M ~ 10 -8 Mq yr _1 for solar-mass 
stars . The relationship has been observed for oyer 200 stars spanning the mass range 0.01 Mq—3 M q 
(e.g., iNatta et alJbood : iMuzerolle et alJbood : Isicilia-Aguilar et al-lbool iMohantv et~aD b()05bl lah . 
There is a large scatter, with values of M for the same M varying by a factor of 10 2 . This scatter is 
intrinsic and not an observational uncertainty. The samples are pre-main sequence stars in young 
clusters. The cluster environment appears to have little effect on accretion: stars in the young, 
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compact p Oph cluster show very similar accretion rates to stars in older, smaller Taurus and Cha I, 
although none of the regions surveyed are as large and dense as is Orion today. 



Figure [12] compares the observations of iMuzerolle et al.l (|2005l ) in Taurus and Cha I (both 
consisting of small cores analogous to our SMALL cluster) to our calculated M-q. We have fit their 
observational data in the same way as eq. [13] and find 



/MX 2 ' 1 
M B (observed) « 21 X 1(T 9 ( — — — 1 M yr' 1 , 
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or about 5 times higher than our computed rate in the SMALL cluster. The slope and scatter clearly 
match well. 

Standard models of disk accretion d o not explain the M oc M 2 relationship, but many new 
models have been prop osed to address it. iGregorv et al.l (|2006l ) improves the standard layered disk 
model (iGammid I1996T) by a s sumin g a temperature structure that depends on the central star's 
luminosity. lHartmann et al.l (|2006l ) assumes a complex magnetic field geometry that introduces a 
mass dependence in the accretion rate. Although these models both move in the correct direction 



neithe r provides a sufficiently steep mass dependence in the exponent. lAlexander fe Armitage 
( 20061 ) develop a model which reproduces the M 2 dependence, but requires certain disk initial 
conditions: namely, that the disk mass scale as the stellar mass (or faste r) , and that the disk radius 
be highest for low-mass stars. Alternatively, IClarke fc Pringld ( 20061 ) argue that the M oc M 2 
dependence may be an observational bias caused by the masking of low accretion rates by the 
bright continua of high- mass stars because the detectability of accretion signatures scales as M 2 . 



Padoan et al.1 (|2005l ) proposed that the observed luminosity could be due to BH accretion from 
the molecular cloud environment. Indeed, they noted that the functional form of BH accretion 
reproduces the M oc M 2 relationship, and the scatter in the observations can be explained by 
variations in stellar speed and gas density. They performed 3D adaptive mesh refinement (AMR) 
calculations of a single star traveling through gas filaments in a Taurus-like region, and found that 
the BHA rate could roughly match the observed accretion rates. 

Our calculated rates are ~5 times lower than those observed, although within the errorbars 
of both, many points overlap. The rough similarity we find between the data and the BHA rates 
is interesting, but it must be interpreted with caution because the observations detect disk-to-star 
accretion, while the BH rates we calculate are for cloud-to-disk accretion. These are two distinct 
physical processes. The two accretions may be ultimately related, buffered through transport 
within the disk. The disk viscous timescale can be estimated by t v ~ r 2 jv for disk radius r and 
viscosity v. Typical vis cous times range from 10 5 yr (r = 10 AU,a = 10~ 2 ) or less, up to 10 7 yr 



(r = 100 AU, a = 10" 3 ; iHollenbach et all 120001 ). If BHA deposits mass at distances of 10-50 AU, 



then this material would be accreted onto the star and cause an observable signature during the 
timescale of star formation. If mass is deposited much further out, or onto disks of low a, then 
BHA' effects would not be observed while the star was still in the cluster. (Or alternatively, the 
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effects of BHA may be manifested as accretion onto the star long after a star has left its molecular 
cloud environment.) 

We raise the possibility that accretion onto the disk may trigger a larger accretion onto the 
star, thus linking these two rates. In particular, because the accreted material can have essentially 
arbitrary angular momentum and can be highly episodic, it may have a destabilizing effect on 
local regions of the disk, causing rapid inward migration and accretion. Almost certainly, multiple 
process contribute to the observed accretion rate and disk dispersal (such as those mentioned 
above), and BHA is only one of these. More detailed modeling of the accretion flow must be 
performed to constrain the dynamics and timescales, and determine to what extent the observed 
and BH accretion rates may be related. 



5.2. Observational signatures of BH accretion 

Bondi-Hoyle accretion will impact the disk surface with speeds of order the higher of the stellar 
velocity or keplerian velocity, typically several km s _1 . Over most of the surface area of a disk, 
this impact velocity will result in shock speeds that are small and post-shock layers that are cool. 
Ongoing accretion from the inner disk onto the central star manifests itself as UV excess, veiling 
at visual wavelengths, and Ha emission; on the other hand, BH accretion onto the outer disk is 
likely to produce a signature visible only at far-infrared and sub-millimeter wavelengths due to 
the accretion-shock velocity. Shock speeds of a few km s _1 will result in post-shock temperatures 
of only a few times 10 2 to 10 3 K. This post-shock layer may be detectable in the pure rotational 
transitions of H2 (such as the 28, 17, and 12 /im lines), in vibrationally excited transitions of species 
such as CO and CS, or high- lying pure rotational states of common molecules such as CO and OH. 
The fine-structure lines of common ions and atoms such as O I, O III, N II, and C II may also be 
excited. High angular resolution observations should reveal that these tracers originate from an 
extended region of the disk comparable to the gravitational accretion radius. Thus, the emission 
produced by BHA should be distinguishable from processes related to the accretion of matter from 
the inner disk onto the star, or from shocks produced by outflows. 



6. Discussion 

BH accretion can provide a substantial amount of material to young disks at a critical time. 
The implications of this process have not been considered by standard models for the formation 
and evolution of disks, planets, or stars. In this section, we discuss several possible implications of 
BH accretion. 
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6.1. Bondi-Hoyle accretion vs. other environmental effects 

Historically, models for disk evolution (including the proto-Solar nebula) have not considered 
the influence of the environment. In the last decade, models have begun to study photo-evaporation 
and close stellar encounters as two potentially significant effects of the environment. First, photo- 
evaporation is a disk dispersal mechanism that rem oves gas and small dust grains by a UV-heating 



mechanism at rate s of 10 -10 M^ yr (e.g.. iThroop &, Ballyi 120051 : iMatsuyama et alJ 12003 ; 



Throop et al.ll200ll ; I Johnstone et al.lll998l ). In large clusters such as Orion the UV source can be 
external O/B stars (which remove the disk outer edges), while in small clusters a star's own UV 
luminosity is the dominant removal source (from the inner edge outward). Second, close approaches 
between stars can strip protoplanetary disks to roughly 1/3 of their miss-distance. These encounters 
rarely affect the inner disk, but can be important for the largest disks in the densest cl usters: typical 



close - approach distances are ~ 500 AU after several Myr in an Orion-like cluster (jAdams et al 



2006) 



Unlike these two processes which remove disks, BHA causes an initial mass increase to the 
disk. BHA may have a wider diversity of manifestations than the other environmental effects. For 
instance, the photo-evaporation rate for far-ultraviolet (FUV) flux depends on disk size, but is es- 
sentially independent of stellar distance. Thus, even though the intercepted flux ch anges quickly for 



stars , they experience minimal change in photo-evaporation through their orbits (jScallv &: Clarke 



200ll ). The BHA rate, however, varies by 10 times or more through a star's orbit, by a few times 
more based on the star's initial position and velocity, and by another factor of (M/Mq) 2 for different 
stars. BHA also occurs in clusters of all sizes, unlike photo-evaporation and close approaches, which 
both occur primarily in the largest and densest of clusters. Initially BHA causes a mass gain, but 
depending on the flow dynamics, the mismatch in angular momentum between the disk and ISM 
may trigger a disk instability, resulting in net mass loss from the disk onto the star. The total mass 
accreted by BHA is comparable to (or slightly less than) the mass loss due to photo-evaporation; 
both of these are in general much larger than mass loss from close encounters. 

Depending on the environment, photo-evaporation and BHA may both occur together, but at 
alternate times. For instance, young stars near an OB association may routinely pass between an 
H II region (where photo-evaporation removes the disk), and a cool molecular cloud region outside 
the Stromgren sphere (where BHA can restore the disk). This 'batter-fry-batter' scenario would 
be possible today in a region like Orion. The Orion Molecular Cloud (OMC) is a region of dense 
non-ionized gas (n = 10 5 cm -3 ) lying just 0.1-0.2 pc behind the Trapezium cluster and the Orion 
nebula. Stars orbiting through the Trapezium region may regularly pass through the OMC, causing 
alternate episodes of photo-evaporation and disk accretion. 

Because of the diversity of effects of tail-end BHA, the full consequences on the formation of 
stars, planets, and disks cannot yet be judged. Most important is an understanding of how the 
inflow interacts with the existing disk, including where and when mass and angular momentum are 
deposited onto the disk. This is likely to be a complex function of (at least) the disk size and orien- 
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tation, and the speed, density, and angular momentum of the local medium. For instance, accretion 
in a filamentary molecular cloud may have a different episodic nature than a much smoother cloud 
or one with lower turbulence. The rotation of a cloud may make a difference, as could its particular 
shape and temperature structure. 



6.2. Stellar mass growth 

As mentioned in § [TJ the role of BHA in the origin of the IMF is the subject of some debate. 

While we have argued that the accretion is predominantly onto the disk and not the star, much of the 

accreted mass will eventually accrete onto the star, because in the absence of photo-evaporation, 

inward accretion is the dominant clearing mechanism for material inward of 100 AU on 10 7 yr 

i II I to J 

timescales (jHollenbach et al.l 2000). An upper limit to the stellar mass increase is simply the mass 



AMb accreted onto the disk. 

Our small model is of a similar mas s and diameter cluster as that in the 'competitive accretion' 
model studied by lBate &: Bonnelll (|2005l ). They found that stars grew via BHA from small initial 
condensations to normal IMF in a few 10 5 yr. In contrast, our results find that BHA is generally 
insufficient to change the mass of most stars: only a handful of solar-mass stars increased their 
mass by as much as 10% during 4 Myr. However, the two simulations are really not comparable, 
because we have si mplified dramati c ally t he interaction between gas and the youngest stars. The 
SPH simulations of lBate fc Bonnelll (|2005l ) model the gas dynamics and filamentary structure to a 
high resolution. Most of the accretion occurs in these filaments, where fragmentation begins and 
the relative velocity between stars and gas is very low. The star stays in these dense filaments 
during most of its formation. Once it is formed, it is often ejected by interaction with another 
young star, or the filament itself disperses. We model a later phase, where stars have dynamically 
decoupled from structure in the local gas. We have taken more conservative estimates for the gas 
temperature and the gas-stellar relative velocity, and we allow for gravitational expansion as the 
gas is removed. We believe that our values are very safe lower limits for accretion rate, but we 
find no reason to believe that the stellar mass function changes substantially during the 'tail-end' 
accretion phase. 

The competitive accretion model predicts very rapid star formation - at most a few 10 5 yr 
- which must be followed soon by cluster dispersal in order to stop accretion from continuing. If 
stars are shown to form on such short timescales, then the tail-end accretion amounts computed 
here would nece ssarily be affected; however, measurements for cluster lifetimes do currently favor 
long timescales (ILada fc Ladall2003l ). 
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6.3. Stellar accretion and metallicity enhancement 

For solar- mass stars and smaller, BH accretion is not enough to substantially increase the 
stellar mass. The accreted mass can affect the stellar metallicity, however, if the composition of 
the accreted material differs from that of the star. This effect can be particularly important for 
stars which are not fully convective. While the lowest mass stars remain fully convective for tens 
of Myr, stars above a few solar masses will develop radiative cores surrounded by a thin convective 
shell within a few Myr of birth. For example, at an age of about 2 Myr, a 2 .5 M star's convective 
envelope contains only the outer 10% of its mass (ILaughlin fc Adamslll997l ). Accretion of high- (or 
low-) metallicity 'polluted' material can lead to enhancement (or depletion) of the star's metallicity. 

Metallicity variations between a star and the accreting material could occur for several reasons. 
First, the cloud itself can vary in metallicity both spatially and temporally. A single supernova 
can significantly alter the metallicity of the surrounding ISM impacted by its ejecta. For example, 
depending on the explosion energy, the demise of a star with an initial mass of 23 Mq can produce 
3 to 4 M m of oxygen, 05 M m of carbon, 0.5 M of silicon, 0.001 M & of 26 A1, and 10 -4 M Q 
of 60 Fe ( Young fc Fryeii 120071 ). If this material were to be mixed with 1O 4 M0 or less of ISM, 
the abundances of the most common elements can be more than doubled and the abundances of 
short-lived radioactive species can be enhanced by orders of magnitude. However, X-ray studies of 
supernova remnants such as Cass-A show that the supernova ejecta is highly inhomogeneous and 
contains clumps with highly non-uniform composition. When such clumps impact the surrounding 
molecular cloud, they can induce even larger local variations in composition. The passage of a 
young, low-mass star surrounded by a disk through such a chemically enhanced region of a cloud can 
result in the accumulation of a veneer of freshly synthesized products of supernova nucleosynthesis. 
This accretion can introduce abnormally high concentrations of common elements such as O and 
C, and short-lived radioactive species such as 26 A1 and 60 Fe. Such elements can then contaminate 
the disk, the stellar photosphere, or both. 

Contamination by superno yae ejecta has been sugge sted to explain metallici t y var iations of 
stars in Orion. Observations by ICunha et al.1 (|2000l . Il998l ) and ICunha k, Lambertl |l994j) of 29 B, 
F, and G stars found abundance variations up to 4 times between stars within the same Orion 
subgroup. These metallicity variations were spatially correlated, lending credence to the idea that 
recent supernovae have contaminated distinct regions of the cluster. The abundance variations were 
seen in O and Si (which are produced by massive stars and type II supernovae), but not in Fe, C, 
and N (which are produced to a much lesser degree, and thus would not be expected to contaminate 
the stellar atmospheres). The observations alone cannot say whether the metallicity variation was 
incorporated into the stars after their formation, or whether it is indicative of heterogeneous mixing 
in the par ent cloud, bu t both scenarios argue for recent injection of large amounts of metals into 
the cloud (|Harpeiill996l ). 



Spatial variation of a molecular cloud's ga s:dust ratio can also lead to metallicity enhance- 
ments in accreted material. IPadoan et al.1 ( 20061 ) used measurements of visual extinction and 13 CO 
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abundance in Taurus to infer variations in the gas:dust ratio of up to 2-5 times within individual 
clouds, on size scales ~ Rb- They also showed that turbulent concentration processes within the 
ISM could lead to such clustering, and that the expected enhancement factors rose quickly across 
smaller distances, to ~ 100 times on AU scales. Accretion over long scales could reduce this con- 
centration, but spatial heterogeneities may persist in the disk for long tim escales if radial mixing 



of small grains within the disk can be inhibited (e.g., iJohansen et al.1 120071 ). 



Supernova contamination requires consideration of the evolution of OB associations. The most 
massive stars can explode within 3 Myr of their birth; the more common massive stars have main- 
sequence lifetimes of 10-30 Myr. A variety of observations have shown that O B associatio ns are 



assembled on a timescale of about 10-20 Myr. Orion presents a nearby example (|Ballvll200ll ). The 
oldest subgroup, Ori-OBla, has an age of 12 to 15 Myr, while the lb and lc subgroups have ages 
of order 3 to 8 Myr. The OBld subgroup is actively forming today and includes the Orion Nebula, 
NGC 2024, and a number of other smaller clusters. These subgroups form a well-ordered sequence 
in both space and time, indicating sequential star formation. 

An important feature of most OB associations, including Orion, is that all of the subgroups, and 
the surviving portions of associated molecular clouds, are located in the interiors of their respective 
supershells and bubbles. The superbubbles are created by the combined effects of ionization, stellar 
winds, and supernova explosions, and are therefore heavily contaminated with supernova ejecta. 
Gamma-ray observations have provided direct evidence for t he presence of large amounts of liy e 



26 Al and 60 Fe in the interior of the Orion / Eridanus feature ( Dullemond et al.l 120071 ; iDiehll 120021 ) . 
The Orion / Eridanus Loops form a 100 by 300 pc Ha emission-line feature that marks the current 
ionization front of the Orion OB association. The energetics of the HI shell that lies outside the 
Ha-emitting rim indicates that at least a dozen supernovae have exploded in Orion over the last 5 
to 10 Myr. All of the surviving molecular clouds and active sites of star formation are embedded 
within the interior of the supershell, and are therefore exposed to the contents of the superbubble. 

Two sites in Orion provide tangible examples of regions that could become SN-contaminated. 
First, the Orion Nebula itself is located directly behind the Ori-OBlc sub-group that contains the 
3-5 Mya clusters NGC1980 and NGC1981. Although there is no evidence for a recent supernova in 
Ori-OBlc, there are massive stars in this complex. When they do explode, they will be in excellent 
position to contaminate the Orion Nebula region and the proximal portions of the L1641 cloud. 
Second, the <r-Ori sub-group of Ori-OBlc is located about 30' or 5 pc west of the southern part 
of the Orion B molecular cloud. Although there is an HI and Ha bubble surrounding this region, 
there is no concrete evidence for a recent SN. However, as seen from the center of the u-Ori cluster, 
the cloud subtends at least 1 steradian, and when an explosion happens, it will pollute the western 
side of the Orion A cloud. Young stars forming from the polluted material, or older stars with 
disks that undergo BHA, will become contaminated with supernova debris. 

Tail-end accretion has a testable consequence for large clusters where metallicity enhancement 
from supernovae is expected to occur after several Myr. Veneers of metal-rich ejecta accreted onto 
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stars in the cluster will have the largest effect on the oldest moderate mass stars which have the 
smallest convective zones. In this scenario, we expect a correlation between increasing metallicity, 
increasing mass, and decreasing stellar age. 



6.4. General effects on Solar System formation 

The planets in the Solar System are believed to have begun formation soon afte r formation 
of the disk. Models build the gas giants within 5-10 Myr (JLissauer fe Stevenson! 120071 ) . consistent 
with the stron g observational req uirement that by 10 Myr, virtually all stars have lost their ob- 

For the terres trial planets, observations indicate that grains 

and 



servable disks ( Calvet et al. 



200' 



begin to grow rapidly in even the youngest disks (jDominik et al. 



2007 



Throop et al. 



2001 



models consistently pr edict runaway growth to planetesimals of size ~ 0.1 Mgarth within 1-10 Myr 
(JNagasawa et alj|2007l . and references therein). Because this is in the same time that tail-end BHA 
is occurring, accretion can affect not only the initial conditions of planet formation, but also the 
physics during the process. 

Tail-end accretion provides a fundamentally new input to the process because it delivers mate- 
rial to the disk well after the disk has formed. The amount of material accreted onto and processed 
by the disk might exceed (by a few times) the original disk mass, or its instantaneous mass. (Model 
disk masses for our Solar System are usually greater than 1 MMSN, but rarely exceed 10 MMSN 
because such large disks become rapidly unstable.) 

This new material may replenish - either locally or globally - gas and dust in the disk that have 
already been used. It may have dynamical effects on the existing disk, because of the mismatch 
in angular momentum. And, because accreted material comes from elsewhere in the ISM, any 
composition heterogeneities in the nebula may be inherited by the disk. 

One 'cartoon' scenario provides an interesting taste of the possibilities. Consider the case of a 
pre-planetary system that begins to form in the Trapezium region, but before O/B stars turn on. 
Its initial disk is several MMSN, and in its first 1-2 Myr it forms a giant planet core. When the 
O/B stars are born, photo-evaporation begins and rapidly removes its remaining gas disk, before 
the planet can build its gas envelope. Subsequently, the star's orbit takes it through the nearby 
OMC region, where BH accretion rapidly replenishes the disk: the 1 Mj (0.1 MMSN) atmosphere 
for the planet can accrete onto the disk in 10 5 yr. This gas can then accrete onto the planet, 
where it is safe against photo-evaporation even after it re-enter s the H II reg ion. Thus, a birth 
environment thought to be destructive to giant planet formation (JThroopll2000l ) may prove instead 
to have regenerative powers. The architecture and compositional history of such planets may lie 
far outside those considered by standard models. 
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6.5. Compositional heterogeneities in the Solar System 

Although it has historically been assumed that the protosolar nebula was well-mixed, recent 
discoveries have shown that this is decidedly not the case. Variations in isotopic and bulk chemical 
abundances in the Solar System strongly argue that it formed from distinct reservoirs of material. 

Isotopic measurements of various species provide extremely precise determination of the chem- 
ical components that formed the young solar system. Isotopic differences between the Earth, Mars, 
and numerous diverse classes of asteroids has provided consistent evidence that multiple reservoirs 
of material were present during formation of the terrestrial bodies. To date, isotopic anomalies 
have been measured for barium, chromium, sul f ur, titanium, zirconiu m, molybdenum, and oxygen 



(JTrinquier et al. 



2007 



Ranen fe Jacobsenl 120061 ; iDauphas et alj|2002l . and references therein). The 



anomalies are small (<3000-10,000 ppm) but definite. They are in general not explai nable by chem- 
ical processes, with the possible exception of oxygen (JYurimoto fc Kuramotd 12004 ) . Importantly, 
the variations are not simply a random scatter, but are correlated. The correlations can be easily 
matched with species produced in nucleosynthetic reactions. Thus, these measurements are usually 
interpreted as evidence for heterogeneous mixing between multiple nucleosynthetic-produced reser- 
voirs during the Solar System's formation (c.f. previous references). While some of these species are 
s-process, and thus produced during the lifetime of massive stars, others are r-process and probably 
produced only in supernovae. 

The observed heterogeneity could be indicative of material accreted onto the disk (e.g., from a 
site spatially distinct from the formation location), or it could be present in the original disk. How- 
ever, studies of short-lived radionuclides such as 60 Fe suggest that at least some of the heterogeneity 
was in fact accreted after the solar system formed. 



6.6. Short-lived radionuclides in the Solar System 



Bizzarro et al.l (J20071 ) compared radioactive decay products present in the Solar System's pri- 



mordial bodies. Iron meteorites and pallasites were formed at < 1 Myr after the start of the solar 
system (where t = is the formation time of the calcium-aluminum-rich inclusions, CAIs). These 
bodies show no evidence for short-lived 60 Fe (£1/2 = 1-5 Myr) incorporated at their formation, based 
on present-day Ni isotopic abundances. However, bodies that formed slightly later at t = 2 — 3 Myr, 
such as the Earth, Mars, and carbonaceous chondrites, consistently show signatures implying they 
were formed in the presence of 60 Fe. Because the only natural sources of 60 Fe are supe r novae , 
and 60 Fe can be incorporated into an exposed disk (JQuellette et al.1 120071 ). iBizzarro et al.1 (J2007I ) 
conclude that this species was not present at initial formation of the protoplanetary disk. Rather, 
it was absorbed later from the environment, presumably a dense cluster with stars of M > 8 M & 
which can produce supernovae. Further, the authors find that 26 A1, another short-lived species 
(ti/2 = 0.7 Myr) primarily produced by massive stars, was present from the Solar System's origins, 
and not added later. This suggests that multiple supernovae in the Solar System's first several Myr 
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contributed material to the young protoplanetary disk. 

The simplest way for supernova ejecta to be incorporated into a disk is by direct implantation. 
Ouellette et al.l (J2007I ) noted that gas-phase supernova ejecta would be deflected by the bow shock 
that forms as the blast-wave over-runs an obstacle such as a circumstellar disk. They suggested 
that most supernova debris is injected in solid form by grains that move ballistically through the 
shock, an d implant in the proto-pla netary disks after the proto-Sun formed (the 'aerogel' model). 
However, IWilliams &: Gaidod ([20071 ) modeled the evolution of stars in clusters of various sizes and 
found the aerogel model to be inefficient: it requires a very massive supernova progenitor, very close 
to the nascent Solar System. They estimated the probability for any disk to acquire the observed 
amounts of short-lived species to be less than a few percent. 

The BH accretion scenario proposed here may provide a viable acquisition route for the short- 
lived radionuclides. First, in an OB association such as Orion, grains of supernova debris can be 
implanted into surviving portions of molecular clouds. YSOs that form from such material initially 
would be enriched in supernova ejecta. Additionally, as pre-existing young stars pass through these 
polluted environments, their disks can accrete a veneer of contaminated gas and dust, adding to an 
otherwise uncontaminated system. Large solids (primitive meteorite particles) that formed before 
passage through the contaminated zone would be devoid of these products, while those that formed 
later would be enriched by them, thus explaini ng the apparent heterogeneity of the decay products 



of short lived species (e.o. lBizzarro et al.ll2007l ). This model has the efficiency advantages of using 



the large cross-section of the entire ISM as a 'dam' for such radionuclides, and using the larger 
'fishing net' radius of Rb vs. th e disk physical radius to c apture material. Thus, the difficulty in 
accreting such nuclides posed by IWilliams fc Gaidosl (J2007I ) may be addressed. 



6.7. Damping of terrestrial planet eccentricities. 

In order to collisionally grow, planetesimals in the terrestrial planet zone require high eccentric- 
ities for orbital crossing. However, the planets today have much lower eccentricities (<2% for Earth 
and Venus) which cannot easily be explained without a damping mechanism. It has been suggested 
that the eccentricities could be damped by a tenuous remnant gas disk of lifetime 10 6 — 10 yr 



(JAgnor fc Ward! 120021 : iKominami fc Idall2002l ). For instance, the latter paper studied a gas disk 
of surface density ~ 10~ 3 — 10 -4 that of the MMSN, which exponentially decayed on several Myr 
timescales. The initial gas density was sufficient to maintain circular orbits. As the disk dispersed, 
orbital interactions entered a chaotic regime which allowed for planetary mergers and growth. The 
remaining gas slowly damped the eccentricities of these widely spaced larger bodies until a small 
number of planets were left on circular orbits. The position and size of the planets depended on 
the amount of gas and its decay profile. 

Tail-end accretion fits into this scenario naturally, and can allow for orbital damping on 
timescales other than the disk's initial dissipation time. The lower limit for the planet formation 



25 



timescale in our sol ar system corr esponds roughly to the 10 Myr upper limit for cluster lifetimes. 
The simulations by iRuffertl ( 1999 ) showed that BHA easily produced disks with 10~~ 3 MMSN den- 



sity, more than sufficient for orbital damping. In a cluster, where the accretion disks are transient 
and episodic, this may provide multiple stages and/or speeds of damping, as a star passes through 
the cluster core multiple times. As terrestrial planets begin to be detected around other systems, 
their eccentricities may give insight into their birth environment as well. 



7. Conclusions 

'Tail-end' Bondi-Hoyle accretion of molecular gas onto young star-disk systems moving through 
the gravitational potential of three realistic model star clusters and their molecular clouds are 
investigated using N-body simulations. We make the following findings: 

• Bondi-Hoyle accretion can occur as stars orbit through dense cloud cores located near their 
birth environments. Gas is accreted toward the star, but is intercepted by the disk before it 
hits the star. Accretion can occur for the several Myr between formation of the first YSOs 
until dispersal of the molecular cloud. 

• For solar-mass stars in Taurus-sized clusters (N=30), the accretion rate is ~ 1O -8 M0 yr^ 1 
(one MMSN per Myr); the rate scales with stellar mass as Mb oc M 2 - 1±ai . 

• Accretion onto stars up to 10 Mq is robust against stellar winds, outflows, radiation pressure, 
and turbulence under most conditions. 

• For a fixed cluster diameter, the accretion rate is weakly dependent on the cluster density. 
Increasing N by a factor of 10 2 decreases accretion rate by only a factor of 3. The lower gas 
density in small clusters is offset by the ease of accretion at their lower stellar velocities. 

• Accretion is an episodic process. The accretion rates during brief plunges through the cluster 
core are up to 10 times higher than the lowest accretion rates for the same star. Intrinsic 
variations in gas density and velocity contribute further to the episodic nature of accretion. 

• The scatter in time-averaged mean accretion rates for stars of the same mass in the same 
cluster is a factor of ~ 7; the scatter in instantaneous accretion rates for the same stars is a 
factor of ~ 50. The scatter is due to the spread in positions and velocities of individual stars 
and by the variations in the density of ambient medium. 

• The mean accretion rate is a few times higher than the median rate, because most accretion 
occurs where stars are moving quickly and least likely to be observed. 

• Bondi-Hoyle accretion onto disks cannot be solely responsible for the M oc M 2 accretion 
observed onto young stars, but these two accretion processes may be linked by mass transport 
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through the disk. The BH accretion rate has similar slope and scatter to stellar accretion 
rates, but is ~ 5 times lower than observed accretion rates. 

• Accretion onto the disk is a low-energy process and difficult to detect directly, but may be 
visible in mid-IR rotational and/or vibrational transitions of H2, CO, CS, and OH. 

• Variations in chemical abundances in the gas cloud can be caused by pollution from massive 
stars. Accretion of this gas can result in a disk with substantially different metallicity than 
its star. Isotopic and bulk chemical heterogeneities in meteorites and terrestrial material may 
be caused by tail-end BH accretion of recent nucleosynthetic products produced by massive 
stars and/or supernovae. Accretion may also explain observed large metallicity variations in 
the atmospheres of young stars. 

• While we estimate the amount of mass deposited during tail-end accretion, the detailed 
dynamics of accretion onto an existing disk are not well understood, and additional modeling 
is needed to understand the deposition of mass and angular momentum onto the disk. 
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Parameter 


Small 


Medium 


Large 


N 


30 


300 


3000 


Number of runs 


200 


20 


2 


M gas 


30 M 


300 M Q 


3000 M Q 


T 


5K 


15 K 


40 K 


Peak stellar density (N pc~ 3 ) 


to 4 


10 5 


10 6 


Peak gas density (n cm- 3 ) 


to 4 


10 5 


10 6 


Stellar mass range 


O.f - 1.2 M 


0.1 - 18 M Q 


0.1 - 27 M Q 


Plummer radius vq 


0.22 pc 


Star formation efficiency (SFE) 


33% 


Mean stellar mass 


0.5 M 


Gas dispersal delay time 


2 Myr 


End of simulation 


4 Myr 



Table 1: Initial conditions for the three cluster models. 



Result, main runs 


Small 


Medium 


Large 


AM B (M ) 


0-038 {£) 


0-018 (^) 


0.008 (^) 


Mean M B (M yr" 1 ) 


9-3 x 10-9 (4-)" U 


4-3 x 10-9 (^ 


2-1 x 10-9 (M.) 


Median M B (M Q yr" 1 ) 


4-7 x lO-^ (^ 


1-6 x 10-9 [M.) 


0.56 x 10-9 (^) Zi 


Result, T=25 K test case 


Small (25 K) 


Medium (25 K) 


Large (25 K) 


AM B (M ) 


0-027 (M-f 


o-oi7 (%y- L 


0-008 (M-) 


Mean M B (M yr" 1 ) 


6-5 x 10-9 (^ 


3-9 x 10-9 (M.) 2 - 1 


2 - lxl0 " 9 te)" 


Median M B (M yr" 1 ) 


3-6 x 10-9 (fiSf* 


1-5 >< 10 " 9 fe)"'" 


0.56x10-9(4)"' 



Table 2: Results for the three cluster models. The top table contains our final results; the lower 
table investigates the effect of cluster temperature on accretion. 
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Bondi-Hoyle Accretion Disk 



Circumstellar Disk 



Fig. 1. — Cartoon of Bondi-Hoyle accretion. The star is moving relative to the ISM gas. Material is 
focused toward the star. Most material accretes behind the star, where opposing streamlines meet 
and cancel their velocities. Accretion is onto a large, temporary accretion disk of size ~ 0.2 R&, 
which deposits material onto the star's circumstellar disk of size < 0.1 Rr. 
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Fig. 2. — Comparison of observed cores to Plu mmer spheres. Th e lower curve plots the column 
density of a pre-stellar c ondensation in L168 9 (jAndre et al.ll2000l ). while the upper curves show 
dense cores in Cassiopeia (IBeuther et al,ll2002l ). All are fit well by the column density of a Plummer 
sphere. 
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Fig. 3. — Overhead views of the large cluster simulation, at start (top) and end (bottom). The 
circular aperture indicates the tq radius of the Plummer sphere; initially, 40% of the mass of stars 
and gas is within this radius. The two magnified views (right) show just the inner region. Larger 
dots indicate more massive stars. 



36 



0.6 

0.4 

0.2 

0.0 

-0.2 

-0.4 

-0.6 



CL 

>- 



Q. 

>- 



Small cluster, 0.997 M 



■_ Small cluster, 0.815 M 








.,...,... 




0.6 


Medium cluster, 


1.19 M 


\ 


0.4 






\ 


0.2 






\ 


0.0 






: 


0.2 






\ 


0.4 






\ 


0.6 










Medium cluster, 


1.03 M 




:..... V ^- 




. . .■ 



— ' — I — ■ — ' — ' — I — ' — ' — ' — I — ' — ' — ■ — I — ' — ' — ■ — I — ' — ' — ' — I — ■ — ■ — ' — I — 

Large cluster, 1.17 M Q 




Small cluster, 0.528 M Q 




Medium cluster, 0.814 M 




■ i ■ ■ ■ i ' ' ' i ' ■ ■ i ■ ■ '\i 
Large cluster, 0.971\M o 




-0.6-0.4-0.2 0.0 0.2 0.4 0.6 
X [pc] 



-0.6-0.4-0.2 0.0 0.2 0.4 0.6 
X [pc] 



-0.6-0.4-0.2 0.0 0.2 0.4 0.6 
X [pc] 



Fig. 4. — Stellar position for nine sample stars, taken from the three simulations. Top row: three 
stars from small cluster. Center: medium cluster. Bottom: large cluster. Stellar masses are 
indicated on the plots, and range from 0.5 — 1.4 M & . Each star's final position is shown with a filled 
circle. Because the orbits are projected into two dimensions, eccentricities appear exaggerated. 
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Fig. 5. — Stellar distance from core. Stars are the same as identified in Figure [H Stars travel 
on highly eccentric orbits. Orbits in the small cluster are the slowest, because of the cluster's 
low potential. As the gas disperses, orbital periods and distances increase, and some stars become 
unbound. 
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Fig. 6. — Stellar velocity and BH accretion rate M-q. Dotted lines: stellar velocity. Velocity 
increases as stars approach the cluster core, and is fastest for the large cluster. Black lines: BH 
accretion rate, based on the 'virial' velocity method (see text). Accretion rate usually correlates 
with velocity, because regions near the core have a higher gas density. Accretion slows after gas 
loss begins at 2 Myr, due to reduced gas density and larger stellar orbits. Small high-frequency 
perturbations are due to close encounters with other stars. Stars are the same as identified in 
Figured! Grey lines: BH accretion rate, based on the 'nearest neighbor' velocity method. This 
is more accurate, but the correlation with stellar velocity is not as obvious. 
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Fig. 7. — Total mass accreted per star (i.e., integral of the solid line in Figured!) Accretion is often 
episodic, with most of the mass being deposited during short periods. Accretion in each 'stairstep' 
can be many Jupiter masses (lMj ~ O.OOIMq). Stars are the same as identified in FigurelU Black 
lines: Virial method. Black lines: Nearest-neighbor method. 
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Fig. 8. — Total mass growth for stars in the three model clusters. Each point corresponds to one 
star in the simulation, and measures its total mass growth by t = 4 Myr. The solid lines are fits to 
the model results, which find AMb oc M 2 ' l±0A across a range of 100 or more in mass. The grey 
grid-lines indicate growth of one MMSN (0.01 M & ) for a solar-mass star. Mass growth is highest 
for the SMALL cluster, although the difference between the three clusters is small compared to the 
scatter within the cluster. Radiation pressure (ignored here) cuts off accretion for the few stars 
above 10 M . 
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Fig. 9. — Instantaneous ('median') values of Mb- Each point represents a single star, at a single 
timestep, and the vertical lines correspond to the varying accretion rates of individual stars through 
their orbits. Solid lines are numerical fits, which find M ~ M 2 ' 1±0,1 . 
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Fig. 10. — Spread in Mb, for solar-mass stars in the large cluster. Histogram shows the range of 
values of M-q for the entire simulation, with one data point per star, per timestep. The la spread 
in accretion rates is a factor of 50. 
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Fig. 11. — Same as Figure [101 but plot shows total mass accreted by solar-mass stars, la spread is 
a factor of 7. 
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Fig. 12. — Comparison between ob servation of acc r etion onto young stars, and calculated BH 
accretion rates. Data shown is from iMuzerolle et al.l ()2005l ) (MLB05); solid line is our fit to their 
data. The dotted line is the result from BH accretion in the small model. The shaded area shows 
the la scatter of M-q in the model. The scatter from the two distributions clearly overlap, although 
the center of the model is several times lower than the data. 



